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^£ ' Abstract 

On 

Radiative QCD corrections significantly influence the theoretical pre- 



dictions for the decay rates of the Z and the Higgs boson. The status 
of the QCD calculations to the hadronic Z width is reviewed. The role 
of mass corrections from bottom quark final states is emphasized. An 
estimate of the theoretical uncertainties is given. New results for quar- 
^ ■ tic mass terms of order 0(a 2 ) are presented. The impact of secondary 

radiation of bottom quarks on the determination of T(Z — > bb) is 
discussed. Second order QCD corrections to the partial decay rate 
T(H — > bb) are also presented in this talk. A recent result for the 
flavour singlet contribution to this quantity is presented. It includes 
quark mass effects and completes the otherwise massless calculations 
of order 0(a 2 ). 
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1 Introduction 

Since experiments at the e + e~ storage ring LEP started data taking a few years 
ago, more than 7 million events have been collected at the Z resonance [0]. Besides 
the electroweak sector of the Standard Model, LEP provides an ideal laboratory 
for the investigation of strong interactions. Due to their purely leptonic initial 
state events are very clean from both theoretical and experimental point of view 
and represent the "golden" place for testing QCD. From cross section measure- 
ments as well as the analysis of event topologies the strong coupling constant can 
be extracted. Other observables which are measurable with very high precision 
are the (partial) Z decay rates into hadrons Ph a d/r e = 20.763 ±0.049 and bottom 
quarks r^/Thad = 0.2200 ± 0.0027 |l[. The program of experimentation at LEP 
is still not completed. The prospect of additional 60 pb _1 per experiment means 
for example that the relative uncertainty of the partial decay rate into b quarks 
ATb/Tb may become even smaller than one percent and an experimental error 
for a s of 0.002 may be achieved ||. 

Also at lower energies significant improvements can be expected in the ac- 
curacy of cross section measurements. The energy region of around 10 GeV just 
below the BB threshold will be covered with high statistics at future B meson 
factories. The cross section between the charm and bottom thresholds can be 
measured at the BEPC storage ring. These measurements could provide a precise 
value for a s and -- even more important --a beautiful proof of the running of 
the strong coupling constant. 

In view of this experimental situation theoretical predictions for the various 
observables with comparable or even better accuracy become mandatory and 
higher order radiative corrections are required. Significant improvements in these 
calculations have been achieved in the last years and will be reviewed in this talk. 
The discussion includes singlet as well as nonsinglet corrections to the vector 
and the axial vector correlator. The corresponding ratios R v and R A enter both 
the formulae for the e + e~ annihilation cross section and the hadronic Z boson 
width. Estimates of the theoretical uncertainties, based on a variation of the 
renormalization scale, will be reviewed in Section 2. Particular emphasis is put on 
mass corrections from bottom quark final states. The quadratic mass corrections 
are of relevance for Z decays and at lower energies. Quartic mass terms, however, 
are of particular importance for the low energy region, especially below and above 
the bb threshold. Corrections of order m 2 /s and new results of order a^m 4 || will 
be reviewed in Section 3. The precise measurement of the rate T(Z — > bb) allows 
for a determination of the mass of the top through its impact on the effective Zbb 
vertex. Secondary bb radiation as well as the assignment of singlet contributions 
are covered in Section 4. 



Figure 1: Renormalization scale dependence of the nonsinglet (a) and the axial 
singlet (b) massless QCD corrections. (a s (Mf ) = 0.12) 

Despite the remarkable confirmation of the Standard Model by the precision 
tests at LEP and SLC, experimental evidence for the Higgs boson is still missing. 
Future colliders like LHC and NLC may provide enough energy for the production 
of the Higgs particle and open the experimental door to the electroweak symmetry 
breaking sector of the Standard Model. In Section 5 QCD corrections to the Higgs 
decay rate are discussed. With T(H — > bb) being the dominant decay channel for 
intermediate Higgs masses, second order singlet corrections to this quantity are 
calculated in the heavy top mass limit. Leading quark mass terms contribute to 
the same order as nonsinglet graphs in the massless limit and are numerically 
sensitive on the ratio M^/ml- 

2 The Z Boson Width 

2.1 Results for Massless Quarks 

Higher order QCD corrections to e + e~ annihilation into hadrons were first calcu- 
lated for the electromagnetic case in the the approximation of massless quarks. 
Considering the annihilation process through the Z boson, numerous new fea- 
tures and subtleties become relevant at the present level of precision. An impor- 
tant distinction, namely "nonsinglet" versus "singlet" diagrams, originates from 
two classes of diagrams with intrinsically different topology and resulting charge 
structure. The first class of diagrams consists of nonsinglet contributions with one 
fermion loop coupled to the external currents. All these amplitudes are propor- 
tional to the square of the quark charge. QCD corrections corresponding to these 
diagrams contribute a correction factor which is independent from the current 
under consideration. Singlet contributions arise from a second class of diagrams 
where two currents are coupled to two different fermion loops and hence can be 
cut into two parts by cutting gluon lines only. They cannot be assigned to the 
contribution from one individual quark species. In the axial vector and the vector 
case the first contribution of this type arises in order 0(a 2 s ) and 0{a^) respec- 
tively. Each of them gives rise to a charge structure different from the nonsinglet 
terms. 

The nonsinglet terms have been calculated in Q and |J to second and third 
order in a s respectively. As shown in Figure la the scale dependence of the re- 



Figure 2: Mass Corrections to R v (left) and R A (right). The left bar represents 
the result in the on- shell-scheme, the right one is obtained in the MS-scheme. 

suit (evaluated in the MS-scheme and for a renormaliation scale fi 2 between s/4 
and 4s) is reduced considerably through the inclusion of higher order terms. The 
dependence of the result due to the renormalization scale is an indicator for the 
influence of yet uncalculated higher order QCD corrections. The remaining varia- 
tion of Rns of about ±2- 10 -4 for a s (M§) = 0.12 translates into an uncertainty in 
a s of ±0.7 • 10 -3 . This error estimate is in itself strongly dependent on the actual 
value of a s . For a s (M|) = 0.135 for example, it is amplified to ±0.001. These 
numbers are quite comparable to the effect of an (a s /7r) 4 term with a coefficient 
around 100 as advocated in ||. 

Similar observations apply for the singlet contributions to the axial current 
correlator (see Figure lb). The order 0(a 2 ) term [0] still exhibits a sizable scale 
dependence. Inclusion of the 0(a^) term ||, however, leads to fairly stable an- 
swer. It should be emphasized that the contribution of this singlet term, which is 
effectively present in the axial bb rate only, is larger than the total contribution 
of the a 3 term. 



3 Mass Corrections 



3.1 Quadratic Mass Corrections 

The calculation of higher order QCD corrections with massive quarks may be 
simplified for small masses by expanding in m? / s <C 1. The operator product 
expansion of current correlators, including subleading terms, provides the theo- 
retical framework. It has been developed in (I7|, [18] and applied to the present 
problem in [IB, H|. The expansion is simultaneously performed in a s and the 
quark mass: 
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The calculation is conveniently performed in the MS-scheme with the running 
mass, which has the remarkable property that no logarithms of the large ratio 



s/m appear in q for arbitrary orders in a s . (This does not hold true for the 
coefficients dj.) For the vector induced rate the first coefficient Cq vanishes and 



the corrections c( , c% , c\ have been calculated in fl(|, fTTH pl and [|13|] respectively. 
For nf = 5 one obtains 



4 1} = 12^^ {l + 8.74^ + 45.15 (*)') (2) 
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In Figure 2 it is demonstrated that the calculation is fairly stable and higher 
order corrections are small in the MS-scheme, in contrast to the large changes in 
the on-shell scheme. 

Similar considerations apply for the axial current induced two point function. 
In particular one can again demonstrate the absence of large logarithms. The ex- 
pansion starts in this case with Cq already. Nonsinglet terms have been evaluated 
in [[TJj to order 0(a^), the leading singlet terms of order 0(a 2 s ) were calculated 



in 15, 16 : 



R (i) = _a< a s 




(3) 



As shown in Figure 2 the calculation in the on-shell scheme exhibits significant 
changes with the inclusion of coefficients with large logarithms. In the MS-scheme 
the expansion is remarkably stable. The size of the corrections is comparable to 
the anticipated experimental precision. 

3.2 Quartic Mass Corrections 

The second order calculation of quartic mass corrections presented below is based 
on 0. The operator product expansion included power law suppressed terms up 
to operators of dimension four induced by nonvanishing quark masses. Renormal- 
isation group arguments similar to those employed already in JC1 [TJ|] allowed to 



deduce the a^m 4 terms. The calculation was performed for vector and axial vec- 
tor current nonsinglet correlators. The first one is of course relevant for electron 
positron annihilation into heavy quarks at arbitrary energies, the second one for 
Z decays into b quarks and for top production at a future linear collider. Below 
only the results for Ry will be presented. 

QCD corrections to the vector current correlator in order a s and for arbitrary 



m / s were derived in flQfl. These are conventionally expressed in terms of the pole 



2 Note that an erratum for the £(3) term turned out to be wrong and the originally published 
result proved to be correct (see 113]). 



Figure 3: Comparison between the complete 0(a s ) correction function (solid line) 
and approximations of increasing order (dashed lines) in m 2 for vector induced 
rates. 



mass denoted by m in the following. It is straightforward to expand these results 
in m 2 /s and one obtains 
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The approximations to the correction function for the vector current correlator 
(including successively higher orders and without the factor a s /ir) are compared 
to the full result in Fig|| As can be seen in FigfJ, for high energies, say for 
2mfe/ \fs below 0.3, an excellent approximation is provided by the constant plus 
the m 2 term. In the region of 2m/\/s above 0.3 the m 4 term becomes increasingly 
important. The inclusion of this term improves the agreement significantly and 
leads to an excellent approximation even up to 2m/\/s ~ 0.7 or 0.8. For the 
narrow region between 0.6 and 0.8 the agreement is further improved through 
the m 6 term. The logarithms accompanying the m 4 terms can also be absorbed 
through a redefinition of m in terms of the MS mass ||iT| P0f at scale s and one 
obtains (m = m(s)) 

Rv = 1-6^-8^ (5) 
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This resummation is possible for the second and fourth powers of m in first order 
a s and in fact for m 2 corrections to Ry and Ra in all orders of a s . However, 
logarithmic terms persist in the m 4 corrections, starting from 0(a 2 ). 

Motivated by the fact that the first few terms in the m 2 /s expansion provide 
already an excellent approximation to the complete answer in order a s , the three 
loop corrections have been calculated. The calculation is based on the operator 
product expansion of the T-product of two vector currents J M = vrf^d. Here u 
and d are just two generically different quarks with masses m u and m^. Quarks 
which are not coupled to the external current will influence the result in order 



Figure 4: Contributions to R v from m 4 terms including successively higher orders 
in a s (order a Q s / a\/ a 2 corresponding to dotted/ dashed/ solid lines) as functions 
of 2m po i e / v /s. 



Figure 5: Predictions for R including successively higher orders in m . 

a 2 through their coupling to the gluon field. The result may be immediately 
transformed to the case of the electromagnetic current of a heavy, say, t (or b ) 
quark. 

The asymptotic behaviour of the transverse part of this (operator valued) 
function for Q 2 = —q 2 — > oo is given by an OPE of the following form. (Different 
powers of Q 2 may be studied separately and only operators of dimension 4 are 
displayed.) 

i [T(J^x)J+(0))e^dx = ^J2(W»-M 2 ) C n (Q 2 ,fi 2 ,a s )O n + ... (6) 

Only the gauge invariant operators G 2 , rriiqjqj and a polynomial of fourth order 
in the masses contributes to physical matrix elements. The coefficient functions 
were calculated in |T7j , the vacuum expectation values of the relevant operators 



in [£lj [18], Q. Employing renormalization group arguments the vacuum expec- 
tation value of J2 n C n O n is under control up to terms of order a s as far as the 
constant terms are concerned and even up to a 2 for the logarithmic terms pro- 
portional to hiQ 2 /fi 2 . Only these logarithmic terms contribute to the absorptive 
part. Hence one arrives at the full answer for a 2 m 4 /s 2 corrections. Internal quark 
loops contribute in this order, giving rise to the terms proportional to J2 m2 an d 
Y^ m i below. 

The result reads (below we set for brevity the MS normalization scale /j, = ^/s 
and m u (s) = m^s) = m; a = a s /ii) 

Tf) / | j 
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Note that the sum over i includes also the quark coupled to the external 
current and with mass denoted by m. Hence in the case with one heavy quark 



u of mass m (d = u) one should set J2i =i = 1 an d J2i =% = 1- In the opposite 
case when one considers the correlator of light (massless) quarks the heavy quark 
appears only through its coupling to gluons. There one finds: 
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The Z decay rate is hardly affected by the m 4 contributions. The lowest order 
term evaluated with m = 2.6 GeV leads to a relative suppression (enhancement) 
of about 5 x 10~ 6 for the vector (axial vector) current induced Z — > bb rate. 
Terms of increasing order in a s become successively smaller. It is worth noting, 
however, that the corresponding series, evaluated in the onshell scheme, leads 
to terms which are larger by about one order of magnitude and of oscillatory 
signs. The m 4 correction to T(Z — > qq) which starts in order a 2 is evidently even 
smaller. From these considerations it is clear that m 4 corrections to the Z decay 
rate are well under control -- despite the still missing singlet piece -- and that 
they can be neglected for all practical purposes. 

The situation is different in the low energy region, say several GeV above the 
charm or the bottom threshold. For definiteness the second case will be considered 
and for simplicity all other masses will be put to zero The contributions to R v 
from m 4 terms are presented in Fig.|] as functions of 2m/\/s in the range from 
0.05 to 1. As input parameters m po i e = 4.7GeV and a s {m 2 z ) = 0.12 have been 
chosen. Corrections of higher orders are added successively. The prediction is 
fairly stable with increasing order in a s as a consequence of the fact that most 
large logarithms were absorbed in the running mass. The relative magnitude of 
the sequence of terms from the m 2 expansion is displayed in Fig.[5[ The curves 
for m° and m 2 are based on corrections up to third order in a s with the m 2 term 
starting at first order. The m 4 curve receives corrections from order zero to two. 

Of course, very close to threshold, say above 0.75 (corresponding to \fs below 
13 GeV) the approximation is expected to break down, as indicated already in 
Fig.|3|. Below the bb threshold, however, one may decouple the bottom quark and 
consider mass corrections from the charmed quark within the same formalism. 

Also Rgg where q denotes a massless quark coupled to the external current 
is affected by virtual or real heavy quark radiation. The m 2 terms have been 
calculated in [[U| and start from order a^: 
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The terms of order a 2 m 4 were given above. Both lead to corrections of (9(10~ 4 ), 
(evaluated at an energy y/s of 10 GeV) and can be neglected for all practical 
purposes. 



4 Secondary bb Production 

The formulae for the QCD corrections to the total rate i\ad have a simple, unam- 
biguous meaning. The theoretical predictions for individual qq channels, however, 
require additional interpretation. In fact, starting from order 0(a 2 s ) it is no longer 
possible to assign all hadronic final states to well specified qq channels in a unique 
manner. For definiteness we shall try to isolate and define T(Z — > bb). The thor- 
ough understanding of QCD and mass corrections is mandatory to disentangle 
weak corrections with a variation of about 1% for m t between 150 and 200 GeV 
from QCD effects. 

1. Non Singlet Contribution 
The vector and axial vector induced rates receive (non singlet) contributions from 
the diagrams, where the heavy quark pair is radiated off a light qq system. The 
rate for this specific contribution to the qqbb final state is given by [p3[] 
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(10) 
Despite the fact that b quarks are present in the four fermion final state, 
the natural prescription is to assign these contributions to the rate into the qq 
channel. Therefore those events with primary light quarks and secondary bottom 
quarks must be subtracted experimentally from the partial rate T b i. This should 
be possible, since their signature is characterized by a large invariant mass of the 
light quark pair and a small invariant mass of the bottom system, which is emitted 
collinear to the light quark momentum. If this subtraction is not performed, the 
bb rate is overestimated by (for m\, = 4.7 GeV and a s = 0.115 . . . 0.18) 

E r NS_ 
1 qqbb 

A = q=uAs ' c ~ 0.005... 0.013 (11) 

Tfefe 

2. Singlet Contribution 
The situation is more complicated for the four fermion final state from the singlet 
contribution which originates from the interference term R S - h i between the qq 
and bb induced amplitudes. It cannot be assigned in an unambiguous way to an 
individual qq partial rate. 



For the vector current induced rate and after phase space integration this term 
vanishes as a consequence of charge conjugation (Furry's theorem). 

For the axial current induced rate, however, this interference term gives a non- 
vanishing contribution which remains finite even in the limit m q —>■ 0. Neglecting 
all masses one obtains 

dS___IdS_ _ i Ip-S __ _ IpS 1 __|_IpS _ _n 1 ko f^A 
^bbbb — 9 bbuu ~ "i" 9 bbdd ~ 2 bbcc ~ "r 9 bbss ~ u,1JO I J 

Among the final states from these singlet diagrams with at least one bb pair, 
only the bbbb term remains after all compensations have been taken into account. 
Numerically it is tiny. Let us compare it with the other singlet contributions 
consisting of bb(g) jets. They can be assigned to r b g in a unique way, although 
they are induced by the bb and the ti currents. These two cuts, in particular the 
leading term with two bb jets, dominate the four fermion final state by a factor 
of 20. Therefore it is suggestive to assign the total singlet contribution to T^. 

5 The Higgs Decay Rate 

The partial width T(H — » bb) is significantly affected by QCD radiative cor- 
rections. First order 0(a s ) corrections were studied in P^fl . Besides the overall 
ml dependence due to the Higgs-bottom coupling, the otherwise massless correc- 
tions in second order 0(a 2 ) were obtained by |[25|| . The resulting expression for 
the decay rate in question reads 
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(12) 



with ATi = f , Ar 2 = 29.14. 

In this talk we present the calculation of an additional contribution to the 
"massless" 0(a 2 s ) corrections, which were neglected in the literature so far. These 
singlet corrections are due to the "light-by-light" type diagram which we calculate 
for nonvanishing quark masses in the heavy top limit. The Yukawa couplings of 
the fermions yield a factor m^nit and each fermion trace produces a factor nib 
and m t respectively. Together with a power supression of 1/m 2 on dimensional 
grounds all mass factors are combined to m 2 , which is of the same order as the 
otherwise massless 0(a 2 ) nonsinglet corrections. 

The absorptive part of the triangle diagram is given by 
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Here also the decay mode into two gluons is included. It is separately known 



26 and should be subtracted in order to arrive at the singlet correction for the 



bottom final state: 
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At last, combining (|12|) and (|14D we get 
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It should be finally stressed that the separation between the decay channels into 
bottom quarks and gluons leads to a logarithmically enhanced In (mf/Mjj) cor- 
rection to the decay width. 
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